Emerging Zero Modes for Graphene in a Periodic Potential 
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We investigate the effect of a periodic potential on the electronic states and conductance of 
graphene. It is demonstrated that for a cosine potential V(x) = Vb cos (Go 2), new zero energy 
states emerge whenever Jo(j 



v F G 



0. The phase of the wavefunctions of these states is shown to 
be related to periodic solutions of the equation of motion of an overdamped particle in a periodic 
potential, subject to a periodic force. Numerical solutions of the Dirac equation confirm the existence 
of these states, and demonstrate the chirality of states in their vicinity. Conductance resonances 
are shown to accompany the emergence of these induced Dirac points. 

PACS numbers: 73.21.-b,73.20.Hb,73.22-f 
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The laboratory realization of graphene [H, 0, [![, a 
two-dimensional honeycomb lattice of carbon atoms, has 
motivated intense theoretical and experimental inves- 
tigations of, among many other things, its transport 
properties [J]. Graphene differs from conventional two- 
dimensional electron systems in that it supports two in- 
equivalent Dirac points, through which the Fermi en- 
ergy passes when the system is nominally undoped. The 
chirality of the electron wavefunctions near the Dirac 
points severely limits electron backscattering, enhancing 
the conductivity of the system The unusual trans- 
port properties associated with the Dirac points may in 
principle be explored and exploited if the edge structure 
can be controlled @, Q , or by application of non- uniform 
potentials, such as in a pn junction d, [l(| EH, EH- 
Periodic potentials may be induced by interaction with 
a substrate [HI, EH EH] or controlled adatom deposition 
[ill. Recently, the existence of periodic ripples in sus- 
pended graphene has been demonstrated [lT}; in a per- 
pendicular electric field this would also induce a periodic 
potential. 

In this paper we discuss the effects of a one-dimensional 
superlattice potential on the transport properties of 
graphene. Previously it has been noted that such a po- 
tential may create a strong anisotropy in the electron 
velocity around the Dirac point [l8(. In fact this be- 
havior is a precursor to the formation of further Dirac 
points in the band-structure of the system. We will show 
that such new zero energy states of the Dirac equation 
are associated with how the phase of the wavefunction 
varies with position along the superlattice direction, and 
is connected with the dynamics of a highly overdamped 
particle subject to periodic potentials both in space and 
time. From an analysis of this problem we demonstrate 
that the emergence of new Dirac points is controlled by 
the parameter Vq/Gq, where Vq is the potential ampli- 
tude (assumed to be a cosine) and Lq — 2t:/Gq is the pe- 
riod. New Dirac points emerge whenever Jo ( ) = 0, 
where Jo is a Bessel function and vf the speed of the 
Dirac fermions in the absence of the potential. The total 



number of Dirac points (associated with a single valley 
and electron spin) is thus 2N + 1, with N the number 
of zeros of Jo(a;), with \x\ < h ^, Go ■ We verify the ex- 
istence of these new zero energy points numerically, and 
demonstrate the chirality of the wavefunctions in their 
vicinity. 

The chiral nature of the wavefunctions around the 
graphene Dirac points in the absence a periodic potential 
has dramatic consequences for its transport properties. 
We find that, when stabilized, the new Dirac points as- 
sociated with the superlattice potential have important 
consequences as well, for transport along the superlattice 
axis. For undoped graphene, strong conductance reso- 
nances appear at the values of Vq/Gq where they first 
appear, as illustrated in Fig. [T] Interestingly, in the 
limit of large width, the conductance scales for most val- 
ues of Vq/Gq as L y /L x , with L y the system width and 
L x its length, indicating the system behaves diffusively. 
This interpretation is supported by the Fano factor, as 
illustrated in Fig. [T^b). At the resonances, however, the 
conductance becomes independent of L x , and the Fano 
factor indicates a more ballistic- like transport. The peri- 
odic appearance of conductance resonances as a function 
of Vq/Gq provides a clear signature of new zero energy 
states as they emerge. 

Zero Energy States - We consider an external poten- 
tial of the form V(x) = Vq cos Gqx, where the period of 
the perturbation is much larger than the graphene lat- 
tice parameter, a, and the amplitude Vq is much smaller 
than the energy bandwidth of the graphene 7r-orbitals. 
In this situation the low energy properties for electrons 
in a single valley and a given spin arc well described by 
the massless Dirac Hamiltonian with a potential, 



H = tiVF(—ia x d x — iuydy) + V{x)X 



(1) 



where a XjV are the Pauli matrices, and X is the identity 
matrix. The wavefunctions which this may multiply have 
two components, $a,s, corresponding to the two trian- 
gular sublattices that make up a honeycomb lattice. 
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FIG. 1: {Color online)(a) Conductivity and (b) Fano fac- 
tor, as function of Vo and for Lo = 50a, for two different 
graphene sample lengths containing 10 and 20 periods of the 
periodic potential VbcosGoa;. Panel (c) shows the conductiv- 
ity as function of Vo / Go as obtained from different superlat- 
tices potentials. This result indicates that the conductivity 
only depends on Vo/Go- Vertical ines indicates the position 
of the zeros of Jo( ~ 2V °„ ) 



To demonstrate the emergence of zero energy states, it 
is convenient to implement a unitary transformation [ 19. 
of the Hamiltonian, H = U\HU u with 



Ui=\ .om .»(.) anda(x) = — 



so that, for &a,b{?) 



— e lkyV c 



V{x')dx', 
(2) 



acting on the transformed wavefunctions Wcj) is 
-idx 



H' = Hvf 



-ik y e ia{xS > 



id x 



(3) 



For 



the cosine potential, 

2Vg \„UGnx 



ia(x) 



YliZ-ooMjd^G- ) e > where J « is the n_th Bessel 

function of the first kind. 



For zero energy states, we search for solutions of H'cj) — 
0. Such solutions have the property 4>a — 4>*b- F°r a sm ~ 
gle valley system this would mean the two components 
are related by time-reversal; this is not quite the case 
here because time reversal in graphene also involves in- 



terchanging valleys. Writing 
single component equation 



6^4 1 e* x , we obtain a 



Or 



y\ + k v e ia - 2i x\<f> A \+i(dx)\<fi A \ = 0. 



(4) 



An analogous manipulation for the zero modes of the 
Bogoliubov-de Gennes equation of states in a vortex of 
a p-wave superconductor plj yields a very simple value 
for the phase x, and the resulting state is a Majorana 
fermion due to the time-reversal relation between com- 
ponents. In our case, we need to find solutions for x that 
are consistent with the symmetries of the problem; when 
they exist the resulting state is that of a real fermion, 
since the two magnitudes 4>a,b are not truly related by 
time-reversal. The equation governing \ is obtained by 
taking the imaginary part of Eq. [H 



k y sin(a - 2x) + d x \ = 0. 



(5) 



The real part yields d x \ 4> a \ + k y cos (a — 2x) 1 4> A \ = , with 
the formal solution 



|c/>a| ocexp| — k y J cos[a(:r') — 2x(x')]riiz;' |. 



(6) 



Since is a Bloch state of the superlattice, it must 
obey the Bloch relation (\>a,b{x + Lq) = e ifca=i ° 4>a,b(x),, 
with k x the crystal momentum. For a zero energy state 
only k x = is possible. We then require (i) \{ x + -^o) = 
x( x ) + 27rm with m an integer, and (ii) J cos[cv(a;) — 
2x(x)] = 0. To sec whether \ can satisfy these relations, 
it is helpful to recast Eq. [5] by writing x = 2x — a, and 
x — > t, so that 



dtX ~ d t a + 2k y sinx = 0. 



(7) 



i,b(x), the Hamiltonian m & ^ 



This is the equation of motion for the position \ of an 
overdamped particle (with unit viscosity), subject to a 
periodic time-dependent force dta and a spatially peri- 
odic force 2k y sinx- Despite the periodicity of the forces 
involved, the generic solution to this equation is not peri- 
odic. However, for certain parameters periodic solutions 
can be found, which correspond to allowed zero energy 
solutions of the Dirac equation in a periodic potential. 
To see this, we solve Eq. [5] perturbatively in k y . Writ- 
0(ky), one finds 



k yX W 



(i) 



dx'sma(x') + C (1) and 

/X pX r>X\ 

dx 1 cosa(x 1 )~ dxi dx 2 sma(x 2 ) + C {2 



where C^ 1 ' 2 ^ are constants of integration. Explicitly per- 
forming the integrations for the above two equations, one 
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finds that condition (i) can be satisfied if 



2G«J L 
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I odd u 



L 



27TTO. 



(8) 



Here Jo and Ji are Bessel functions evaluated at 
2Vo/hvpGo- Since we have employed a small k y expan- 
sion, the only consistent solution is for m = 0. In this 
case Eq. [8] fixes C^, and the resulting x ( an d the as- 
sociated x) 15 periodic. Condition (ii) may then be im- 
plemented to fix the value of k y at which a zero mode 
appears, 



h y 



2 £4, 



'2 odd 



(9) 



Eq. [9] predicts the presence of a zero mode whenever the 
right hand side is positive. This turns out to occur when 
x = is just above the values for which Jq(x) — 

0; we have confirmed numerically that the sign of the 
denominator on the right hand side of Eq. [5] always works 
out such that k y > in this situation. With increasing x, 
the solution moves to larger \k y \ until it diverges where 
the three Bessel function sum vanishes, which is always 
prior the next zero of Jo(x). We note that since our 
approximation is only valid for small k y , Eq. [5] cannot 
accurately predict the location of the zero energy states 
well away from k y = 0. However, since zero energy states 
can only annihilate in pairs [22j, we expect that once 
they emerge from the origin, they should persist. We 
now show that a numerical solution of the Dirac equation 
supports this expectation. 

Numerical Studies - Our expectations about the new 
zero energy states can be directly confirmed by numer- 
ically solving the Dirac equation in a periodic poten- 
tial. To accomplish this we represent the Hamiltonian 
H = Hq + Vo cos Go a; in a plane wave basis and diag- 
onalize the resulting matrix for momenta (k x ,k y ), with 
— Go/ 2 < k x < Go/2. We have checked that our results 
converge with respect to the number of plane wave states 
used. 

Because of the chiral nature of the Dirac quasiparticles 
the effect of the periodic potential is highly anisotropic: 
at the k x = 0, k y — Dirac point the group velocity is un- 
changed in the direction of the superlattice but is strongly 
reduced, even to zero in some cases, in the perpendicular 
direction [H |M H|] . In Fig[2]we plot, for different val- 
ues of Vo , the lowest few energy eigenvalues as a function 
of k y for k x —0 and Lq — 50a. As Vo increases the group 
velocity at the Dirac point decreases to zero [18|, [23j, [2J] , 
and thereafter two zero energy states emerge from k y =0 
as the group velocity of the k y =0 Dirac point becomes 
finite again. These are the new zero energy states dis- 
cussed above; we find that they emerge precisely when 
Jo( ftvG ) = ®- Up 011 further increase of Vq, the group 
velocity along k x at k y —0 becomes zero again and a new 
pair of zero energy states emerge from k y =0, again pre- 
cisely at the next zero of Jq ( nv^-Go ) = ^' This pattern 
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FIG. 2: (Color onZine)Energy bands of graphene in presence 
of a superlattice potential VbcosGoa;, as function of k y with 
k x ~0, for several values of Vo, in units Hvf/ci, and Lq = 50a. 



continues to repeat itself with increasing Vo- Further 
studies for different periodicities confirms the prediction 
that the emergence of these points depends only on the 
ratio Vo/ Go, precisely as predicted in the analysis of the 
previous section. 

These zero energy points in fact represent new Dirac 
points, as can be demonstrated by studying the chirality 
of the wavefunctions in their vicinity. In Fig.?? we plot 
the expectation value of the pseudospin (< a x >, < a y >) 
for the lowest posotive energy eigenvalue as a function of 
crystal momentum. One may see that this vector under- 
goes a 27r rotation for any path enclosing one of the zero 
energy states. The non-vanishing winding of this vector 
for such paths is a clear signal of the Dirac-like nature of 
the spectrum in the vicinity of a zero energy state. 

Conductivity. Using the transfer matrix method we 
have computed the transmission probability through a 
graphene strip of length L x , containing N p periods of 
the superlattice potential. We assume metallic contacts 
connected to the strip may be modeled by heavily doped 
graphene (25| . Boundary conditions are taken to be pe- 
riodic in the transverse direction, leading to transverse 
wavefunctions which can be labeled by a momentum k v ; 
this is justified when the width of the strip, L y , is much 
larger than its length. From the transmission probability 
of each mode, 7\. , we obtain the conductance G and the 
Fano factor (ratio of noise power and mean current) F, 



G = 4YTu F= Efc y fc " (1 ~ r0 (10) 



where the factor 4 accounts for the spin and valley de- 
generacy. The conductivity is related to the conductance 
via geometrical factors, a = G x L y /L x . In what follows 



we work in the limit L y » L x 
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FIG. 3: (Color online) Expectation value of the vector field 
(o~ x , (Jy) of the lowest energy positive subband for a superlat- 
tice of period 50a and Vo = O.lhvF/a. The zero energy Dirac 
points at ^^=0 are highlighted. High energy Dirac points with 
with opposite chirality also appear between the k y —0 Dirac 
point and the superlattice induce Dirac points. 

For pristine graphene, Vq = 0, the conductivity is in- 
dependent of L x and takes the value = 4e 2 /irh, and 
the Fano factor takes on the universal value 1/3. In the 
absence of disorder one might expect the conductivity 
of graphene to be zero or infinity and the the electri- 
cal current to be noiseless. The deviation from these 
naive expectations, consistent with diffusive transport, is 



a unique property of the Dirac points, which has been in- 
terpreted in terms of the spontaneous creation of virtual 
electron- hole pairs [25j . 

FigQ] shows the conductivity and the Fano factor, as 
a function of Vo, for two graphene strips, respectively 
containing 10 and 20 periods of a potential of the form 
Vo cos Go x . For finite values of Vq, apart from some res- 
onances which we discuss momentarily, the system be- 
haves diffusively (F — 1/3) and the conductivity is well 
defined. Interestingly, between the peaks the overall scale 
increases with Vo, showing that the periodic potential 
tends to enhance the conductivity. 

At certain values of Vq one observes peaks in the con- 
ductance, for which the conductivity is not well defined 
and the Fano factor tends to zero. These resonances oc- 
cur precisely whenever new zero energy states emerge 
from the origin in k— space, and represent a direct exper- 
imental signature of their presence. We believe the reso- 
nances occur because the group velocity vanishes when a 
zero energy state emerges, leading to a strong enhance- 
ment of the density of states. A further check that the 
resonances are associated with the zero energy states is 
to see that they depend on the ratio Vq/Gq) Fig. rj|c) 
demonstrates that this is the case not just for the reso- 
nances but for the entire conductance curve. 

In summary, we have shown that a search for zero en- 
ergy states in graphene in a periodic potential shows that 
they may be stabilized for large enough Vb/Go, and that 
their presence may be detected in the conductance of the 
system. The generation of new Dirac points was very 
recently reported in Ref. [26j| . which focuses on different 
properties than the ones described above. 
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